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Phase space path integral is worked out in a riemannian geometry, 
by employing a prescription for the infinitesimal propagator that takes 
riemannian normal coordinates and momenta on an equal footing. The 
operator ordering induced by this prescription leads to the DeWitt cur- 
vature coupling in the Schrodinger equation. 
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1. Introduction 

The research on path integral in curved space began with DeWitt |l|, who 
found that the measure of Pauli-MoreteH, adapted for curvilinear coordinates, 
leads to a coupling with the curvature in the Schrodinger equation (see also 
Ref.p|). Later on Parker M shown that one can take profit of the ambiguous 
integration measure in the configuration space to generate different values for 
the coupling constant in the wave equation. This coupling can be understood as 
a contribution of an ambiguous measure, which can be written in an exponential 
way to dress it up as a potential S,Ti R in the action. The ambiguity in the choice 
of the measure is nothing but the reflection of the different operator orderings 
for the Hamiltonian operator in curved space. 

The ambiguities characterizing path integration methods in curved space, 
can be reexamined in the context of phase space variables. In phase space the 
integration measure seems to be better established: the Liouville measure dq dp 
is privileged because it is invariant under canonical transformations. However, 
as was demonstrated by Kuchafg, ambiguities are still present in phase space. 
Namely, the prescription for skeletonizing the canonical functional action fails 
in being clean, as it is in configuration space. So different prescriptions again 
lead to different curvature coupling in the wave equation (see also Ref.B , and 
Ref.jjj for the handling of a relativistic system). In order to fix the ambigui- 
ties, it has been proposed 8 that the natural skeletonization of the canonical 
functional action should be one managing canonical coordinates and momenta 
on an equal footing. The skeletonized canonical functional action proposed in 
Rcf.[f| is based on those complete solutions of the Hamilton- Jacobi equation 
that accomplish boundary conditions which are symmetric in coordinates and 
momenta. The aim of the present paper is to apply the results of Ref . M to a 
phase space path integral in curved space, and compare the result with DeWitt 
and Parker prescriptions for configuration space path integration. In Section 2 
we review Parker's path integral in configuration space. In Section 3 we present 
the phase space path integral proposed in Ref. M. In Section 4 we apply this 
later prescription to the system studied by DeWitt and Parker. In Section 5 we 
display the conclusions. 



2. Configuration space path integral 

The propagator K(q' t' \ q t) is a two point function that builds those 
solution ty(q' , t') of the Schrodinger equation satisfying a given initial condition 
^o(q) at time t: 

*(q',t')=[dqK(q't'\qt)* (q) (1) 

Thus the propagator itself must be a solution of the Schrodinger equation, ac- 
complishing the boundary condition 

Jim + K(q't'\qt)=S(q',q), (2) 



and the composition rule 

K(q" t" | q' t!) = J K(q" t"\ q t) dq K(q t \ q> t'). (3) 

When curvilinear coordinates are used, one should care about the geometri- 
cal character of the quantities involved in Eq.(pl). In fact, the inner product 

(*i,*2)=[dq*l(q)*2(q) (4) 

must be invariant under general coordinate changes. So the wave functions must 
be regarded as densities of weight 1/2, in order that dq ^l(q) ^2(0) results to 
be an invariant volume. Some people prefer regarding the wave functions as 
invariant; then they add a factor ^fg in the inner product (g is the determinant 
of the metric). The latter choice amounts the replacement \E r — > g 1 ^ 4 ^, and the 
consequent redefinition of the operators. In this paper we will attach ourselves 
to the former choice. Thus, we will regard the propagator as a bidensity of 
weight 1/2. 



The propagator can be formally expressed as a path integral 



K(q t'\ qt) = Vq(t) exp 



\ %(*)] 



(5) 



where the integration is performed on all paths joining the arguments of the 
propagator. This functional integration can be dismantled into infinitesimal 
(short time) propagators K £ (q', q) = K(q' t + e\q t), which can be assembled 
by means of the composition rule (0), so rebuilding the finite time propagator. 
Each infinitesimal piece must be a bidensity of weight 1/2. This skeletoniza- 
tion process give sense to the otherwise ill-defined functional integration (g). 
However the choice for the infinitesimal pieces is not unique. Thus path inte- 
gration is an ambiguous operation: each choice for the skeletonization amounts 
a different quantization of the system, all of them possessing the same classical 
limit. 

The infinitesimal propagator K e (q' , q) is defined by developing the time evo- 
lution of the wave function in the following way: 

Vtf, t) + e JU(g', t) = J d n q K e (q', q) *(q, t) + 0(e 2 ). (6) 

One can recover the wave equation in Eq.dq) by retaining only the terms 
that are linear in e. For this aim one replaces ^(q,t) in the integrand of the 
Eq.(||) by the expansion 

d^ 1 (9\P 



where £ l = q l — q n . 

If the system io governed by the Lagrangian 

L = 2 9ij(<l) <F tf , i = l,...,n. (8) 

then the Schrodinger equation io a second order invariant equation in curved 
space, whose general form is: 



ih^-Mq, t) = -]-h 2 A^ 2 ) *(q, t) + nh 2 R*(q, t) (9) 



where A^ 1 / 2 ) is the Laplacian for densities of weight 1/2: A' 1 / 2 ) = g 1 ! 4 A g~ 4 / 4 , 
R is the curvature escalar of the manifold and /i is a real number. Therefore, 
the terms 0(£ 3 ) in Eq. (7) cannot generate contributions linear in e when they 
are replaced in Eq. (6). In Eq. (9), the coupling constant fi is specified by 
giving a prescription to path integrate the propagator (5) or, equivalently, to 
order the hamiltonian operator. 

Parker's infinitesimal propagator for the system (8) can be expressed in the 
following way || : 

K e (q> | ^^-V ^W'y f gW) -i/2 DeWAq) giq) -i/2V 

(2irih) v ' 





cx]> TSe(q', I q) 
where D{q't' \ q t) is the Van Vleck determinant |1C]] 



(10) 



^YI.^detl-^^M), (11 ) 



and S(q' t' \ q t) is the Hamilton principal function, i.e. the functional action 
S = J L dt evaluated on its stationary (classical) path. 

Parker's propagator (|1C]) satisfies the wave equation 







dt 1 



1*9.1 11/11 1 



ih^-K{q' t'\qt)= - -h 2 A> (V2) + -(l-p) h 2 R(q') K(q' t> \ q t). (12) 



2 6 



The choice p = 1/2 in Eq. ( |l0| ) is the DeWitt's propagator 0], while p = 
corresponds to the Feynman's propagator pi] . 

By composing the infinitesimal pieces ( |10| ) one gets a meaning for the path 
integral (0). Thus the Van Vleck determinants take part in the measure T>q[t) 
and the Hamilton principal functions S (q't' \ qt) are the pieces that skeletonize 
the functional action S[q(t)]: 



AT-1 



%(*)] 



^2 S(q k+1 i fc +i | <?fc tk), 



(13) 



fe=0 



where tk+i — tk = e, and q^ is a shorthand for q l (tk)- 



3. Phase space path integral 

Phase space path integral, 



K(q' t'\ q t) 



Vp(t) Vq(t) exp 



S[q(t),p(t)] 



(14) 



also requires a skeletonized version to make sense. The infinitesimal propagator 
proposed in Ref.flSj is 



K s {q>,q)=K{<ft" = t' + e\q>t')=j-^ 



d n Pa 



h)> 



det 



d 2 J{q"t"\pt) 
dq" j dp a 



1/2 



det- 



d 2 J{pt\q't') 
dp a dq' j 



1/2 



exp 



- {J{q"t"\pt) + J(pt\q't')) 



(15) 



Here J(qt' \pt) and J{pt \qt') are two invariant functions related to the dynam- 
ical system, which should take on an equal footing canonical coordinates and 
momenta. Besides these functions should provide a skeletonized version of the 
canonical functional action S[q(t),p(t)]: 



(JV-2)/2 



S[q(t),p(t)] 



2 / {J(q2k+2t2k+2\P2k+lt2k+l) + J {pik+ltik+l |??2fc*2fc)} ■ 



fc=0 



(16) 

The functions J(qt'\pt) and J(pt\qt') used in Ref.g are two complete 
solutions of the Hamilton- Jacobi equation linked by the relation 



J(qt'\pt) = -J(pt\qt'), 



(17) 



The invariance of J(qt' \pt) is not guaranteed by the invariant Hamiltonian 
in the Hamilton-Jacobi equation, because the boundary condition must also be 
invariant. The boundary condition proposed in Ref. H is 



J(q i t\p a t)=p a a a (q j ), 



(18) 



a a (q J ) in Eg. (Jig) are riemannian normal coordinates for the point P identified 
by the coordinates q l , based at some origin O: 



a a (qi) = s(0,P)u a , 



(19) 



where s(0, P) = J Q ■^g i jdq l dq3 is the (invariant) lenght of the geodesic joining 
O and P, and u a are components of the unitary vector at O that is tangent to 
the geodesic joining O and P, in some basis chosen at O. Riemannian normal 
coordinates have a twofold behavior: like s they are invariant under changes 
of coordinates q l , but transform as components of a vector under changes of 
basis in the target space To- p a in Eqs.fll5|) and (18) are momenta canonically 
conjugated to a a {qi), 

Thus the boundary condition ( p"8| ) is invariant under both general changes of 
the q l, s and change of basis in the target space To- The boundary condition 
(p"q) means that canonically conjugated riemannian normal coordinates and mo- 
menta enter the skeletonization on an equal footing. 

As is shown in Ref.g, J(q' ' t' \pt) and J(p't' \qi), there named Jacobi prin- 
cipal functions, result to be the Legendre transforms of S(a a (q') t' \ u a (q) t) that 
interchange the variables a a and p a at t and t' respectively. Then J(q' t' \pt) 
generates contact transformations 

, dJ(q>> t'\p a t) „ dJ(q'* t'\p a t) 

P * = — d7~> — ' a = — dp~ a — ' (21) 

So the determinants in the infinitesimal propagator (|15| ) are well defined in all 
phase space. The functions J(qt' \pt) and J{pt'\qt) also provides a proper 
skeletonization for the canonical action: 

J(q2k+2 hk+2 \P2k+l hk+l) + J{P2k+l t2k+l\<l2k hk) 
> Pa 2k+l ( a 2k+2 _ a 2k) ~ H\ a 2kiPj 2fe+l,^2fc+l) (hk+2 — hk)- (22) 

Differing from the infinitesimal propagator ( |10|), w hich only holds for systems 
whose Lagrangian is quadratic in the velocities ]l2fl , the infinitesimal propaga- 



tor (15) can be applied to any hamiltonian system. It will provide a unitary 
evolution when t in Eq. (pj) is the mid time t = t' + (e /2) = t" - (e /2). 

Finally we remark that it is by no means evident that the prescription (n5h 
does not depend on the choice of the origin O for defining normal riemannian 
coordinates. This topic will be analyzed in the next Section. 



4. Infinitesimal propagator in the riemannian geometry 

Let be g a b(q) the components of the metric tensor at P in the riemannian 
normal coordinates system (RNCS), 

g ab (q) ^ |£(ff) ££(<,) 9ij (q) g«\q) = ^-(q) ^(q) g ^ {q) (23) 



da a ^' do bV1 ' °" v *' v/ d<t y 'd(t 



g a b is invariant under changes of q l, s, but is a tensor in To- Since u a in Eq.(|l9| 
is unitary (g ab u a u b = 1), one recognizes that 



„ *a ~i> n 2 

g a b o a = s . 



(24) 



Then, 



and 



1 f) <! f) <? 



(25) 



« a (<?)=.9 a "(?)^(<?) 



(26) 



At the lowest order in e = At, the solution of the Hamilton- Jacobi equation 
fulfilling the boundary condition ([18) is 



Mq l I Pa) = Pa<T b (q l ) - \ g a " Pa Pb e + 0( £ 2 ), 



(27) 



d 2 J(e/2)(q\Pa,e) 



dq l dp a 



1/2 



da a 



dqi 



1/2 / h o a ^j 



Both Jacobi principal functions, the one of Eq. (J27|) and its companion of Eq.(O), 
enter the Eq.rtia) to compute the infinitesimal propagator. The result is (see 
Appendix) 



K e {q"t" =t' + e\q't' 



i/2 



(2wih) n/ 



-^ [detg l] (q"T 1,i D 6 (q" | q>) [det grAq'T 1 ' 4 



exp a 



-2 ,t" 



Se (q" I «') - Y / dt (detg ab ) 1/4 A (detg ab )-^ 



(29) 



where the integral on t is evaluated along the classical trajectory q(t) joining 
(q',f) and (q",t"). This infinitesimal propagator leads to the following wave 
equation: 



Q +-2 -1 

ih^(q,t) = -yAl 1 ' 2 ' tffat) + -?i 2 (det 5afc )*A [(det( 5a6 ))-*] *(«,,*). 



The term M Q (?) = % (det g ab (q))i A \(detg ab (q))-i 



(30) 
is invariant under 

changes of the coordinates of point P, where the laplacian operates, and also 
under change of basis in Tq. Except for the factor exp — j? f , dt Mo (q(t)) , 
the result (p9| ) agrees with (|l(j) for p = 1. However this factor depends on the 



choice of the origin O, which is a non admissible dependence for the propagator. 
To be convinced that Mq really depends on O, one can evaluate Mq in RNCS: 



(det gafc) 1 / 4 A (det^)- 1 / 4 



(det g ab 



-1/4 



d_ 



-(det g a b) 



1/2 qfc d(detg ab ) ^ 4 

5 c^ 6 



1 n ab n cd d 9cd ]_ n cd^9cdab n ef^9el _ 1 9 9 n cd 9 9 



i 9 9 da a da b 16 5 dcr Q 



9 



cd 



da b 4 9cr a y da b 
This expression can be compared with the curvature scalar R(q): 

1 dg ab dg de 



R 



2 da c daf 



(31) 



(g ac g de g bf - g ab g de g cf + 2 g dc g be g af + 2 g df g be g ac ) 



d 2 g a b 



(g ac g bd - g ab g cd ) 



da c da dK* » » » , (32) 

Since they do not coincide, the term Mo depends on the choice of the origin 
O to define the p a 's ( in flat space M is null because RNCS is the cartesian 
system where the derivatives of the metric are null) . 



5. Prescription for the origin of RNCS in the infinitesimal propagator 

In order to avoid any undesirable dependence on geometric objects non linked 
to the physical system, the point O in (29) should be dictated by the dynamical 
system itself. So the point O for the infinitesimal propagator (29) must be 
chosen on the classical path joining q"(t") and q'(t'). This prescription will 

give a proper meaning to the factor exp —j-J t , dt Mo(q(t)) in Eq. (f29|). 

The simplest prescription is to choose the point O as the final point in the 



infinitesimal propagator. So Mo{q) in Eq. (30) must be replaced by Mp(q), 
being P the point where the wave function is evaluated. Having into account 
the expansion (||) 



then 



,da a 



d6t( W ] = 9 ~^ P)9 ~^°' ) ( l + « S RlJs > iS 'i + ° ( s ' 



1 



(33) 



limM (q) = —R(q). 



(34) 



Thus the Schrodinger equation(30)turns out to be: 

ih^(q,t) = -ItfA* 1 /*) *(q,t) + ±h 2 R*(q,t) 



(35) 



which exhibits the curvature coupling prescribed by DeWitt. 



6. Conclusions 

To some extent the prescription O = P in the infinitesimal propagator seems 
arbitrary, in the sense that any other point of the classical trajectory joining 
q' and q" is equally natural to play that role. However, the result of Section 5 
-the DeWitt curvature coupling- is not affected by a displacement of point O 
along the classical path. In fact, such a displacement would change M (q(t)) 
to M (q(t)) + 0(£ a ); thus K e (q" \ q') would get a contribution that can be 
expanded in powers of £ a . Because of the rule £ a £ 6 — ► i e h g ab (see Eq. ( p9[) 
in the Appendix), these changes do not modify the leading term in Eq. (6) but 
only contribute to higher order in e. Then, these changes can be ignored at the 
level of the infinitesimal propagator. Therefore the DeWitt curvature coupling 
appears to be the natural one for a quadratic system. 



Appendix 



I Identities 



9ij,k — ~g%h 9 ,k9lj 







a N ■ — Nn N a mn a 


oo 
-oo 


d n x 
(2^ ^ 


-~(x,Ax) + (b,x) 


= (det , 



(36) 

(37) 



(det A)~ 1/2 exp - (b, A^b) (38) 



where A is a n x n matrix and b and x are vectors. 



j Aq^Aq a \..Aq a - exp f^-g^AqiAq^ d n q = (2nhef /2 g~^ 2 (itie) 



{g aia2 ...g a2m ~ ia2m } + all relevant permutations of ai,...,a.2m (39) 



II Derivation of Wave Equation 



The expression (15) can be integrated with the Jacobi principal function 
(27). The result is: 

K 'U' I *) = ; !, n/2 9 ,i 9"hdet(g: b )rHdet( 9 i ab ))-i(det(g ab ))^ 
(2nhej ' 



i lItj ( d a' M d 9" bd ,,- 



4 aa0 V 5cr' d 9(7 



exp ( ^«* 



(40) 



where g afc ee \ (g" ab + g' ab ) . Let be g = a" 3 - a' j . Then: 



9ab k=0= 9ab 



9g ab i 10Sa6 

o— 



<9ct /c 

^ab 



i«=o= 



2 9cr /c 

1 <V ab 

2 <9cr' c 



(41) 
(42) 
(43) 



d 2 g ab | _ 



1 ^ Pat , ! fdg'amdg 



2 <9cr" c <9cr" d 4 V <9cr" c 9a" d 



~9nb "I" 5am 



g g //mn ^ 

<9er" d <9ct" c 



(44) 



By expanding det(g' ab ) and det(g ab ) around <r" a , 



(det^)-i=(det^)-i-^^H3 e+ , 



1 9 2 (det(g^))- t fefd 
2 dcr" c <9er" d 



(45) 



and using the identities (42-44), one obtains the terms that are relevant in the 
integration of Eq. (6). According to Eq.(B9) these terms are: 



(det g ab )i = (det g'^) i (det <4) i + (det <&) * 



1 dg" am dg% 

16 der" d 9^"" c 



re d + 0(C 3 ) (46) 



/6c FirJIbc 



dab £' 



a fdg' bc dg 



V <9cr /c 9ct' c 7 <9cr" c <9cr 



a 2 3," 



-HcftrrllA 



an tine /-a /-d 

9 s s " 



// "3 "9mn„/>nctatd , ^ff ^9an cacd , Q /V3\ 



5a b 



g CT //d 9o ./, 



acr //d ao ./, 



(47) 



exp ^5ah f£ 6 = exp 



a /-b 



2?ie 



5a h re 



* "5afc to rb tc_ 



4he da'' 



10 



pp. II ■ a II Q limn 

u Hab c a C b C c C d _J_ "Jam u i> „ll tagb^cta 

S S S S T ^77 n „,, n »rf 9n& S ? S ? 



8he da" c da" d s s s 8fte da" c da" d 

1 0<& ass 



Jjte r ^ r ^ re / + j (48) 



32ft 2 e 2 <9ct" c 9. 

These expressions are replaced in Eq. (40) to obtain the infinitesimal prop- 
agator K e (q" | q') entering the integration in Eq.(6). This integration is per- 
formed with the help of the identities ( J3^ ) . Thus the wave equation results to 
be 

lh ^J£A = _l^ A(1 / 2) ^ ( ^ i)+ l a 2 ( (det(5 // 6)) 4A(det( ff ;' 6 ))-T) *(«",*) 

(49) 
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